AD-A199  191 


m 


cnF  fAPv 


/  1  /  y  ^ 


1 


TECHNICAL  REPORT  RD-GC-86-16 
VIBRATION  OF  PRESSURIZED  CYLIIOERS 


Carl  H.  Warren 

Guidance  and  Control  Directorate 
Research,  Development,  and  Engineering  Center 
and 

A.  R.  Barbln 
Auburn  University 
Auburn,  Alabama 


<S> 


JANUARY  1988 


Q 


<&!).$.  air  m  mmu lb 

Redstone  Arsenal,  Alabama  35M-5000 


Approved  for  public  release;  Distribution  vnUmited. 


dtio 

ELECTEI 


OCT  04  1988 


\ 


M  FORM  1061,  1  AUG  SB  PREVIOUS  EDITION  IS  OBSOLETE 


88  10  4  070 


DISPOSITION  INSTRUCTIONS 


DESTROY  TIM  IMPORT  WHEN  IT  M  NO  LONGER  NEEDED.  DO  NOT 
RETURN  IT  TO  THE  ORIGINATOR. 


DISCLAIMER 

THE  FINDINGS  M  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OP  THE  ARMY  POSITION  UNLESS  SO  DESIG¬ 
NATED  BY  OTHER  AUTHORIZED  DOCUMENTS. 


TRADE  NAMES 

USE  OP  THAOS  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT  DOES 
NOT  CONSTITUTE  AN  OFFICIAL  INDORSEMENT  OR  APPROVAL  OF 
THE  USE  OF  SUCH  COMMERCIAL  HARDWARE  OR  SOFTWARE. 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

RD-GC-86-16 


lb  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6«  NAME  OF  PERFORMING  ORGANIZATION 

US  Army  Missile  Command 


6c  ADDRESS  (Cry,  State,  and  ZIP  Code) 

Redstone  Arsenal,  AL  35898-5254 


6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

AMSMI-RD-GC 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


8a  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

72W801 


PROJECT 

NO 


TASK 

WORK  UNIT 

NO 

ACCESSION  NO 

00 

7260 

1 1  TITLE  (Include  Security  Classification) 

Vibration  of  Pressurized  Cylinders 


12  PERSONAL  AUTHOR(S) 

Carl  H.  Warren,  Guidance  and  Control  Directorate,  and  A.  R.  Barbin,  Auburn  University 


13a  TYPE  OF  REPORT  13b  TIME  COVERED  114  DATE  OF  REPORT  (Year.  Month.  Day)  15  PAGE  COUNT 

Final  from  Jan  83  TOJul  86  1988  January  20 


16  SUPPLEMENTARY  NOTATION 


COSATI  CODES 


GROUP  SUB-GROUP 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  study  was  conducted  of  the  effects  of  internal  pressure  and  media  on  the  vibration  of 
small  cylindrical  pressure  vessels.  Utilizing  approximate  expressions  for  elastic  poten¬ 
tial  and  kinetic  energies  of  the  cylinder-fluid  system,  LaGrange's  equation  was  applied  and 
an  expression  developed  for  predicting  the  frequency  of  vibration  as  a  function  of  the 
pressure  and  density  of  the  fluid  for  axial  wave  factors  X  =  nmR/L  <  0.60,  thickness  ratios 
6/R  >  0.03  and  number  of  circumferential  waves  n  ^  2.  End  conditions  are  empirically  in¬ 
cluded  in  the  measured  zero-pressure  frequency .  A  series  of  experiments  on  a  small  cyl¬ 
inder  of  the  type  used  as  pneumatic  power  sources  on  missiles  confirmed  the  predictions 
over  a  range  of  pressures  from  0-34  MPa,  with  helium  and  nitrogen  gases.  A  technique  of 
nondestructive  pressure  measurements  is  proposed  based  on  the  results  of  this  study. 


20  DISTRIBUTION  I AVAILA8ILITY  OF  ABSTRACT 

□  UNCLASSIFIEDIUNLIMITED  B5AME  AS  RPT  □  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Carl  H.  Warren 


21  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


OO  FORM  1473,  84  MAR  83  APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


22b  TELEPHONE  (Include  Area  Code)  I  22c  OFFICE  SYM801 

AV  746- 1 225/ 205-876- 122p  AMSMI-RD-GC 


SECURITY  CLASSIFICATION  OF  This  PAGE 

UNCLASSIFIED 


i/(ii  blank) 


TABLE  OF  CONTENTS 


Page 


I.  INTRODUCTION .  1 

II.  ANALYSIS .  2 

III.  EXPERIMENTS .  6 

IV.  RESULTS -DISCUSS  ION .  8 

V.  SUMMARY .  10 

REFERENCES .  11 

GLOSSARY .  12 

APPENDIX  KINETIC  ENERGY  OF  A  FLUID  IN  A  VIBRATING  CIRCULAR 

CYLINDER . .  A- 1 


Accession  For 

NTIS  GRAil 
DTIC  TAB 
Uriflnnoiui-'ed 
Justification _ 


By - 

Di3trlbnt Ion/ 

Aval l.tl  1 1  tty  Codes 


Avail  ani/or 
Diet  Special 


iii/(iv  blank) 


I .  INTRODUCTION 


The  essential  features  of  free  vibration  of  unpressurized  cylindrical 
shells  have  been  well  established  by  previous  investigations  of  infinitely 
long  [1]  and  finite-length  [2,3]  cylinders.  Cylindrical  shells  can  vibrate 
with  various  combinations  of  axial  and  circumferential  waves.  For  any  given 
wave  pattern,  three  distinct  natural  frequencies  are  possible  depending  on 
the  dominant  amplitude  (radial,  axial,  or  circumferential)  of  the  motion  [1], 

According  to  Arnold  and  Warburton  [1],  investigations  of  cylindrical 
shells  can  proceed  either  directly  to  a  set  of  differential  equations  via 
equilibrium  considerations  of  a  small  element  or  indirectly  by  assuming  vibra¬ 
tion  forms  compatible  with  the  given  end  conditions,  deriving  expressions  for 
strain  and  kinetic  energies,  and  using  the  Lagrange  dynamical  equations  to 
arrive  at  the  differential  equations.  Discrepancies  between  the  results  of 
various  investigations  can  be  traced  to  different  approximations  utilized  in 
relating  strains  to  the  extensional  strains  of  the  middle  surface  and  its  cur¬ 
vature  changes;  for  an  accounting  of  the  relations  used  by  various  authors, 
see  the  review  presented  in  Reference  [3]. 

Previous  investigations  of  cylindrical  shell  vibrations  have  been  limited 
to  unpressurized  vessels.  The  purpose  of  this  study  was  to  assess  the  effect 
of  pressurization  on  the  natural  frequency.  The  analysis  is  based  on  a  sim¬ 
plified  energy  method  using  the  assumption  that  radial  displacements  dominate. 
The  zero-pressure  natural  frequency  is  assumed  to  be  given  either  by  existing 
predictive  techniques  or  measurements. 


II.  ANALYSIS 


Figure  I  illustrates  the  cylinder  considered  in  this  treatment.  The  ends 
of  the  cylinder  z  ■  0  and  z  »  L  are  assumed  to  remain  circular  during  vibra¬ 
tion,  and  the  middle  surface  is  presumed  to  move  in  a  strictly  radial  direc¬ 
tion;  thus,  the  displacements  of  the  middle  surface  are 

u  -  U  cos  (n9)  sin  (miz/L)  (1) 

v  «  w  m  0 

where  n  is  the  number  of  circumferential  waves  and  m  is  the  number  of  axial 
halfwaves;  U  is  periodic  in  time 

U  -  UQ  sin  ait  .  (2) 


8 


Figure  1.  Cylinder. 


The  total  potential  energy  of  the  cylinder-fluid  system  is  the  sum  of  the 
cylinder  bending  and  stretching  energies  and  the  fluid  expansion  energy.  In 
computing  the  energy  of  cylinder  stretching,  it  is  assumed  that  tension  is  due 
to  internal  pressure  only,  and  per  unit  length  is  pR  in  the  circumferential 
direction  and  pR/2  in  the  axial  direction;  fluid  pressure  p  is  assumed  con¬ 
stant.  Upon  deformation  the  length  of  the  middle  surface  increases  by 


2 


2  if 

A40  =■  /  (i  (R+u)2  +  (3u/39)2]1/2  -  R)d9  (3) 

o 

in  the  9-direction.  The  integrand  of  this  expression  can  be  expanded  and, 
upon  retention  of  the  first-order  terms  only  (for  small  displacements  u/R  <C  1 
and  slopes  3u/R39<^  1),  approximated  by 

2  it 

Mq  •  /  Iu  +  2r  C-||)2]de  •  (4) 

o 

Similar  treatment  of  extension  parallel  to  the  z-axis  yields 

2  ft 

Aiz  -  f  £  ( 3u/ 3z)2  dz  (5) 


for  extension  of  a  generator  of  the  cylinder.  Equations  (4)  and  (5)  can  be 
used  to  form  the  epxression  for  stretching  energy 


2  TT  L 

/  / 


p[u  + 


2R 


(-f§)2]Rd  9dz 


2  it  L  ^ 

/  /  *S-  <t>2 

o  o 


d  9dz 


(6) 


The  volume  of  fluid  inclosed  by  the  cylinder  increases  during  deformation 
by  the  amount 

2  tr  L 

/  /  uRd 9dz ,  (7) 

o  o 


thus,  the  fluid  expansion  energy  is  given  by 

2 1  L 

vE  -  -  /  /  puRd 9dz  .  (8) 

o  o 

Equations  (6)  and  (8)  can  be  combined  to  yield  the  sum  of  the  stretching  and 
expanison  energies  as 


Vs  +  *E  -  f  f  Ptlg  (f)2  +f  (-g)J]Medr  (9) 

0  0 

or,  the  displacement  of  equation  (1)  can  be  used 

Vs  +  VE  -  x  Ln2U2p[ 1  +  ( A/n)2] /4  (10) 


where 


X  =  mit  R/L 
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The  potential  energy  of  cylinder  bending  Vg  is  calculated  from  the 
expression 

2  it  L  6/2  , 

Vg  -  /  f  f  — - 27~  [ez2  +  e02  +  2 aeze 0  +  (“J2)  Y0z2]Rd9dzdx  (11) 

o  o  -5/2  2(l-<^) 

where,  following  Love  [5],  it  has  been  assumed  that  normal  stress  in  the 
radial  direction,  orr,  is  small  and  the  shear  strains  yrg  and  yrz  are  zero. 
The  strains  ez  and  eg  are  presumed  to  be  linear  functions  of  x,  with  zero 
midplane  strain;  thus  [5] 

ez  -  -x(32u/3x2)  ■  (mir/L)2  Ux  cos  (n8)  sin(mirz/L)  (12) 

e@  »  -x(  32u/  3 g2 ) /r2  >  (n/R)2  Ux  cos  (n0)  sin(mirz/L) 

where  the  displacement  Equation  (1)  has  been  invoked.  These  strain  expres¬ 
sions  may  be  substituted  into  Equation  (11)  resulting  in  the  expression 

<|)3<^2>2  • 

The  total  kinetic  energy  of  the  cylinder-fluid  system  is  the  sum  of  the 
kinetic  energy  of  the  cylinder 

2  it  L  R+6/2 

T c  *  /  /  /  J  P<-^)2  rdQdzdr  -  £  pRL  (U)*6  (14) 

o  o  R-  6/2 

and  the  kinetic  energy  of  the  fluid  (Equation  (A-7).  Thus, 

T  -  Tc  +  Tf  -  itLRpe(U)2  6/4  (15) 

where  pg  is  an  "equivalent  density"  defined  by 

Pe  =  P  +  Rpf/(n6)  .  (16) 

Equations  (10)  and  (13)  may  be  added  to  yield  the  total  potential  energy 
of  the  cylinder-fluid  system;  thus, 

V  “  ^-CT^r  12R^~  (  x2+n2)2  +  i  Ln2u2Pt1+(  X/n>2l  U?) 

Equations  (17)  and  (15)  may  then  be  utilized  in  Lagrange's  equation 

(3T/3U)  -  3T/3U  ■  -3V/3U  (18) 

to  yield  the  expression 
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1 


R5pe 

~2 


(  X2+T12  )  2 

(l-o2) 


+  |  (X2+n2)}  . 


For  periodic  oscillations,  Equation  (2),  this  expression  reduces  to 


<ii2R6pe  -  - ^ — r—  (^)  (X2+n2)2  +  p(X2+n2) 

12(1—0^)  R 


(19) 


(20) 


At  zero  gage  pressure,  p  ■  0,  pe  =  p,  and  Equation  (20)  reduces,  in  that 
instance,  to 


u^2R 6p - -  (|)3  (  X2+n2)2  (21) 

^  Ud-o2)  R 

where  denotes  the  unpressurized  oscillation  frequency.  Equation  (20)  can 
be  conveniently  nondimensionalized  by  division  by  Equation  (21);  after  some 
manipulation  the  results  are 


+  ^  £  A3 

E  6 


(l-o2)  ' 

(X/n)2+l_ 


(22) 


This  expression  can  be  simplified,  with  the  aid  of  Equation  (16)  to  the  form 


0) 


1 

1  R  Pf 
n  1  T 


£(V 

E  6 


(l-o2)  -,)1/2 

(mxR/nL)2  +  1 


(23) 


The  derivation  of  Equation  (23)  is  based  on  the  assumption  that  the  num¬ 
bers  of  axial  and  circumferential  waves  do  not  depend  on  pressure.  The  zero- 
pressure  frequency  is  considered  to  be  given  either  by  experiment  or  some 
predictive  technique  other  than  Equation  (21).  Equation  (23)  indicates  that 
cylinder  pressurization  tends  to  increase  the  frequency  through  the  increased 
strain  energy  with  a  moderating  influence  of  the  kinetic  energy  of  the  fluid. 


5 


III.  EXPERIMENTS 


A  series  of  oscillation  frequency  measurements  were  made  on  a  small  gas 
storage  bottle,  typical  of  those  used  as  pneumatic  storage  supplies  on  certain 
missiles  [4].  Figure  2  shows  the  bottle  and  the  dimensions  of  the  bottle. 
Qualitative  preliminary  measurements  were  conducted  to  determine  the  numbers 
of  axial  and  circumferential  waves  existing  during  free  vibration;  View  A  of 
Figure  3  shows  the  test  setup  used.  An  oscillator  was  used  to  drive  a  small 
loudspeaker  placed  adjacent  to  the  bottle  and  adjusted  to  a  frequency  which 
produced  a  peak  output  of  an  accelerometer  glued  to  the  bottle.  By  moving  the 
accelerometer  axially  and  circumferentially,  and  by  comparing  its  output  to  a 
second  accelerometer  fixed  to  the  bottle,  the  node  lines  were  determined.  Two 
circumferential  waves  (n  ■  2)  and  one  axial  half-wave  (m  »  1)  were  found  to 
exist  for  pressures  from  atmospheric  up  to  34.5  MPa  (5000  psig)  for  nitrogen 
and  helium  gases.  For  the  bottle  tested,  mirR/L  =  0.576  and  5/R  =  0.0867. 
According  to  Arnold  and  Warburton  [1]  these  ratios  indicate  that  the  minimum 
frequency  should  occur  for  n  =*  2,  as  observed. 


Figure  2.  Gas  storage  bottle. 


Frequency  measurements  were  conducted  as  shown  in  the  test  setup  in  View 
B  of  Figure  3.  Vibration  of  the  bottle  was  initiated  by  sharply  tapping  it 
with  a  small  hammer.  The  output  of  a  microphone  placed  adjacent  to  the  bottle 
was  recorded  on  an  oscillograph.  The  measurements  were  completed  by  manually 
counting  the  cycles  in  a  specified  time,  short  enough  (0.100  s)  to  avoid  sig¬ 
nal  dropout  due  to  damping.  Tests  were  made  for  nitrogen  and  helium  gases 
over  a  ringe  of  pressures  from  0-34.5  MPa  gage.  Pressurization  between  mea¬ 
surements  caused  significant  temperature  changes.  Great  care  was  taken  to 
cool  the  cylinder  down  to  room  temperature  (20  °C)  via  a  water  bath  prior  to 
each  measurement. 

Measurement  uncertainties  are  estimated  as  two  parts  in  two  hundred  for 
frequency  ( ~  50  Hz)  and  0.5  MPa  for  pressure.  The  dimensions  of  the  storage 
bottle  shown  in  Figure  2  have  an  uncertainty  of  0.025  mm. 
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VIEW  B.  Frequency  Measurements 


IV.  RESULTS- DISCUS SION 


The  results  of  the  frequency  measurements  are  presented  in  Figure  4  for 
both  helium  and  nitrogen.  The  zero-pressure  frequency  was  3550  Hz  for  each. 
At  the  highest  pressures  for  which  data  was  taken  (  ~  34  MPa),  the  frequency 
for  helium  Increased  by  approximately  26  percent  over  the  unpressurized  case; 
nitrogen  had  a  corresponding  increase  of  16  percent. 


P.MPa 

Figure  4.  Variation  of  frequency  pressure. 


Also  shown  in  Figure  4  is  the  variation  in  vibration  frequency  with  pres¬ 
sure  as  predicted  by  Equation  (23),  using  the  measured  zero-pressure  frequency 
05.  Helium  and  nitrogen  densities  were  calculated  with  a  Van  der  Waals  prin¬ 
ciple  of  corresponding  states  [6].  Agreement  with  the  measurements  is  excel¬ 
lent,  differing  by  less  than  4  percent  of  the  frequency  increases  experienced 
at  the  higher  pressures  and  well  within  the  estimated  uncertainty  of  the  mea¬ 
surements.  Note  that  the  frequency  predicted  by  Equation  (23)  is  somewhat 
sensitive  to  the  thickness  of  the  bottle  5;  for  an  uncertainty  of  0.025  mm  the 
predicted  frequencies  could  have  a  variation  of  +  30  Hz  at  the  higher  pres¬ 
sures. 


The  agreement  between  the  measurements  and  Equation  (23)  might  seem  some¬ 
what  fortuitous  considering  the  rather  restrictive  simplif ications  made  in  the 
analysis.  However,  the  more  comprehensive  analysis  of  Arnold  and  Warburton 
[1]  indicates  that  for  axial-wave  factors  X  <  0.60,  thickness  ratios  5/R  > 
0.03,  and  n  ^  2,  the  strain  energy  in  an  unpressurized  cylinder  is  almost 
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entirely  composed  of  bending,  and  that  the  lowest  frequency  of  vibration  (of 
three  possible)  is  comprised  mainly  of  radial  motion;  these  assumptions  were 
made  at  the  outset  of  this  analysis,  which  should  not  be  extrapolated  beyond 
the  limits  quoted. 

The  conditions  assumed  in  the  analysis  above  differ  from  the  tests  per¬ 
formed  in  one  important  aspect,  the  end  conditions  imposed  on  the  cylinder. 

The  displacements  of  Equation  (1)  are  representative  of  the  condition  of 
"freely-supported”  ends,  since  the  test  cylinder  had  hemispherical  end  caps,  a 
condition  somewhere  between  freely-supported  and  fixed  ends.  In  their  experi¬ 
ments  on  cylinders  with  solid  and  flanged  ends,  Arnold  and  Warburton  [2]  pro¬ 
posed  a  correlation  based  on  appropriately  modified  axial-wave  factors.  Their 
correlation  yields  an  accurate  zero-pressure  frequency  if  the  hemispheric-end 
capped  cylinder  is  considered  as  solid  with  end  thickness  equal  to  cylinder 
thickness  and  with  length  equal  to  the  cylinder  length  plus  end  caps.  In  this 
manner,  is  estimated  at  3560  Hz,  which  is  remarkably  close  to  the  measured 
value  of  3550  Hz. 

Equation  (23)  indicates  that  for  circumferential  waves  n  >  2,  the  effect 
of  pressurization  on  frequency  should  be  small  except  at  very  high  pressures. 
The  notion  that  internal  fluid  pressure  in  a  container  does  not  affect  the 
frequency  [2],  although  a  good  guide  for  most  applications,  is  not  precisely 
correct  and  could  lead  to  significant  underestimations  of  natural  frequencies 
of  high  pressure  storage  vessels  under  appropriate  conditions. 

The  original  impetus  for  this  study  came  from  a  need  to  nondestructively 
find  the  fluid  pressures  inside  small,  high  pressure,  sealed  containers  such 
as  those  used  on  missiles.  A  quick  and  simple  method  is  needed  to  identify 
those  containers  which  have  lost  their  charge  through  leakage  without  time 
consuming  and  costly  disassembly  of  the  system.  This  study  gives  such  a  meth¬ 
od  for  containers  similar  to  those  shown  in  Figure  2.  If  natural  frequency 
measurements  could  be  made  on  the  containers,  pressures  could  be  easily  esti¬ 
mated  to  an  accuracy  of  *0.7  MPa,  sufficient  to  determine  the  operational 
status  of  the  containers. 
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V.  SUMMARY 


This  study  is  summarized  as  follows: 

1.  Internal  pressurization  of  thin  cylinders  can  significantly  increase 
vibration  frequencies  in  those  circumstances  where  the  strain  energy  in  the 
unpressurized  cylinder  is  mainly  bending.  The  effect  is  diminished  as  the 
number  of  circumferential  waves  is  increased  and  is  enhanced  with  decreasing 
fluid  density. 

2.  Natural  frequencies  of  vibration  of  pressurized  thin  cylinders  with 
6/R  >  0.03  and  X  ■  mirR/L  <  0.60  are  well  correlated  by  Equation  (23)  over  a 
range  of  pressures  from  0.34  MPa,  for  helium  and  nitrogen  gases. 

3.  Internal  pressures  in  thin  cylinders  with  6/R  >  0.03  can  be  determined 
via  frequency  measurements  ("ping  test")  to  an  accuracy  of  approximately  2 
percent  at  pressures  up  to  34  MPa. 
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GLOSSARY 


Te  rm 

er ’  ® 6  >  e  z 
E 

Aie 

Aiz 

L 

m 

a 


P 

r,  0  ,  z 

R 

t 

T 


U,  V  ,  w 

U 

V 


x 

a 

5 

Y r9  >  Y0Z  >  Yzr 

v 

V 

\ 

r(a) 

m 

4> 

P 

Pf 

Pe 

o 


Nomenc lature 

Radial,  circumferential,  and  axial  strain 

Young's  modulus,  2.07  x  105  MPa 

Modified  Bessel  function,  first  Rind,  nc^  order 

Circumferential  extension  of  the  middle  surface 

Axial  extension  of  the  middle  surface 

Length  of  cylinder 

Number  of  axial  half-waves 

Number  of  circumferential  waves 

Pressure 

Radial,  circumferential,  and  axial  coordinates 
Radius  of  middle  surface 
T  ime 

Kinetic  energy;  (  )c  and  (  denote  cylinder 

and  fluid,  respectively 

Radial,  circumferential,  and  axial  displacements  of 
middle  surface 

Amplitude  of  radial  displacement 

Potential  energy;  (  )B,  (  )g,  (  )g  denote 

bending,  stretching,  and  expansion,  respectively 

Radial  distance  from  middle  surface 
General  parameter 
Cylinder  thickness 

Shear  strains  in  r-6,  0-z,  and  z-r  planes,  respec- 
t ively 

Distance  normal  to  surface 
Denotes  a  unit  normal 
Axial  wave  factor,  muR/L 
Gamma  function 

Frequency;  (  )Q  refers  to  zero-pressure  case 
Velocity  potential  of  fluid  in  cylinder 
Density  of  cylinder 
Fluid  density 

Equivalent  density  defined  by  Eqn.  (16) 

Poisson's  ratio,  0.286 


Superscripts 

(*)  Denotes  differentiation  with  respect  to  time 

(  )'  Denotes  ordinary  derivative 
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APPENDIX 


KINETIC  ENERGY  OF  FLUID  IN  A  VIBRATING  CIRCULAR  CYLINDER 

The  flow  field  interior  to  the  cylinder  is  assumed  to  be  irrotational  and 
to  behave  incompressibly ,  since  only  small  pressure  variations  are  anticipated 
with  the  small  amplitude  vibrations  considered  here.  The  kinetic  energy  of 
the  fluid  is  given  in  terms  of  the  velocity  potential  $  as 


L  2  if 


Tf  *  J  T?  dA  “  T  pf  f  {  *  T$  rd9dz 


(A-l) 


With  the  approximation  that  the  unit  normal  to  the  surface  is  essentially 
radial,  for  small  displacements,  3$/3v  =  34>/  3r  and,  for  thin  cylinders  with 
6  <$  R,  Equation  (A-l)  can  be  approximated  by 


i  L  2* 

Tf=rpf//  ♦■gwed. 


(A-2) 


The  integrand  Equation  (A-2)  is  evaluated  at  r  »  R,  the  middle  surface  of  the 
cylinder. 

The  function 


-UL  ,  ,  v  ,  ,*n  (mirr/L) 

4.  -  —  sin  (mirz/L)C08(n 9) 


(A-3) 


satisfies  Laplace's  equation  in  the  cylindrical  region  and  the  boundary  con¬ 
dition 


-jjp  (R,  9,z)  »  Usin(nnrz/L)cos(n0) 


(A-4) 


at  the  surface  r  *  R,  in  agreement  with  the  displacement  given  in  Equation 
(1).  Further,  the  net  volume  flow  rate  past  the  surfaces  z  *  0  and  z  a  L, 


2it  a* 

S  f  lwrd9dr 

0  0  32 


(A-5) 


is  zero;  thus,  the  velocity  potential  of  Equation  (A-3)  represents  the  veloc¬ 
ity  field  within  a  cylinder  capped  by  fixed  ends  (of  any  shape). 

Substituting  Equation  (A-3)  into  Equation  (A-l),  results  in  expression 


JT  (U)2L3  In(m*R/L) 
^  4  miT  I^(mitR/L) 


(A-6) 


A-l 


Equation  (A-6)  can  be  simplified  with  the  recurrence  relation 


In(o)  -  In_i(a)  -  £  In(a) 


(A-7) 


and,  for  small  arguments,  the  expression 


In(  °) 


(  ot/2  )n 

'  r(n+l) 


(A-8) 


to  the  form. 


In(mirR/L)  1 

I^(mirR/L)  =  2L  T  (n+1)  nL 

mirR  f  (n)  ni® 


mirR 

nL 


(A-9) 


Equation  (A-6)  can  thus  be  written  as 

Tf  -  ir  pfRL2(U)2/4n  (A-10) 

It  is  interesting  to  note  that  Equation  (A-10)  is  identical  to  the  results 
obtained  by  summing  the  kinetic  energies  of  small  axial  slices  of  the  cylin¬ 
drical  region,  ignoring  completely  the  z-component  of  velocity  in  each. 
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